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Abstract A lemma of Micchelli's, concerning radial polynomials and weighted sums 
of point evaluations, is shown to hold for arbitrary linear functionals, as is Schaback's more 
recent extension of this lemma and Schaback's result concerning interpolation by radial 
polynomials. Schaback's interpolant is explored. 

In his most-cited paper, [M], Micchelli supplies the following interesting auxiliary 
lemma (his Lemma 3.1). 

(1) Lemma. IfYl7=i (^iPi^i) = ^^r all p G n<fc(lR'^), then 

n n 

1=1 j=i 

where equality holds if and only if 

n 
i=l 

Here, (xi, . . . , Xn) is a sequence in IR^, and || ■ || is the Euclidean norm in M'^. Further, 

with 

0" : IR'^ IR : X h-^ := x(l)"(^) ■ ■ ■ 
a convenient if nonstandard notation for the power function, 

n<fc := n<fc(lR'^) := span(()" : a G |a| < k) 

is the collection of all polynomials in d real variables of (total) degree < k, and 

n<fc :=n<fc(iR^) 

those of degree < k. The collection of all polynomials in d real variables will be denoted 
here, correspondingly, by 11 = n(lR'^). 

Micchelli follows the lemma by the following 

(2) Remark. Applying the Lemma inductively shows that the conditions 

n 

^c,p(xO=0, pGn<fc(lR'^), 

1=1 

and 

n n 

(-i)'EE^*^^-^(ll^^-^jH'') = 0' 9en<fc(iR), 

i=i j=i 



are equivalent. 

The essence of the proof is, perhaps, the observation (imphcit in MiccheUi's proof) 



that 



\^-yf' = {\\xf-2j2<i)y{i) + \\yf)' 

i 



(3) ^ a\(5\c\ 

a+\0\+c=k ^ 



with 



= E (-2)' E ^^'".M^KMi/), 



and with a,b,c nonnegative integers, /? e ZZ^, and 

Each summand in the final sum of (3) is thus the product of a (homogeneous) polynomial 
in X of degree 2a + b and a (homogeneous) polynomial in y of degree 2c + b. Hence, if 

i.e., A is a linear functional on 11 that vanishes on n<fc, then the tensor product of A with 
itself, i.e., the linear map 

A ® A : n (8) n ^ IR : 0" 0^ ^ AQ" A()^, 

annihilates all the summands with 2a + b < k or 2c + b < k. As to any other summands, 
they must have 2a + b > k and 2c + b > k, hence 

2k = 2{a + b + c) = 2a + b + 2c + b>2k, 

therefore 

2a + b = k = 2c + b. 

These are the summands in which the polynomial in x equals the polynomial in y and, 
moreover, k — b is even. Thus, altogether (and in faulty but understandable notation), 

I.] ^ 

{X^X)\\x-yr = J2 ^ \((k-b)/2V)^ E (Mfe-^)A/3)V/^!- 

k-b even |/3|=6 

In particular, (A ® A)((— l)'^||x — y\\'^'^) > with equality if and only if A vanishes on the 
sequence 

(4) {P(k-b}/2,t3 :\l3\=b;be k-b even). 

But since each P(k-b)/2,i3 here is homogeneous of degree k while, in particular, k — b is even 
when b = k hence each with \(3\ — k appears in (4), this last condition is equivalent to 
having A vanish on span(()'^ : \P\ = k, P e ^+), hence on Il<k- 

Altogether, this proves the following generalization of MiccheUi's Lemma and Remark. 



(5) Proposition. If X ± n<fc, then 

{-l)''{X^X)\\x-yf^ >0 

with equaUty iff X _L n<fe. 

(6) Corollary. A _L n<fe if and only if 

{X® X)\\x-yf'' = 0, 0<r<k. 

Note the following useful 

(7) Corollary. The bihnear form 

(, )fe : n' n' : (A, //) ^ (-1)'=(A ® f,)\\x - yf" 
on the algebraic dual II' of 11 is an inner product on any algebraic complement of 

±n<fe :={iieU' : ii± U<k} 

in -Ln<A;. 

Schaback [S] reiterates Micchelli's results and extends them as follows (Lemmata 8 
and 9 of [S], though only for linear functionals A that are linear combinations of point 
evaluations) . 

(8) Lemma. If A ± U<k, then, for all k < 21, 

x^X\\x--f^ 

has degree < 2i — k. 

Conversely, if, for some k < 2£, 

x^X\\x--f^ 

has degree <2£ — k, then X ± Il<k- 

To be sure, the first assertion follows directly from the basic identity (3) since, by that 
identity, application of such A to — kills all summands with 2c + \P\ < k, leaving 
only those with 2i - 2a - \/3\ =2c+ \/3\ > k, i.e., with 2i - k > 2a + \/3\. 

For the second assertion, simply "apply the same idea as in the proof of Micchelli's 
lemma", to quote [S]. Arguing perhaps differently, rewrite (3) in terms of polynomial 
degrees in x to get 

j=0 2a+b=j |/3|=6 '^^ 
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Now, (—2)^ = {—iy2'' since j — b here is always even. Also (using a + b + c = £), 
Pa,p = II • ||^^-'"^Vc,/3- Therefore, 

2£ 

j=0 2a+b=j |/3|=6 '^^ 

Hence, if now x i— > A||a; — -Ip^ is of degree <2i — k, then each of the sums 



2a+b=j 1/3 1 =6 



must be zero, hence so must be the value of A on each such sum, i.e., 

2a+b=j |/3|=6 '^^ 

This implies that 

Apc,^ = 0, \/3\ = b, 2£ - k <2a + b, £ = a + b + c, < a, 6, c, 
hence, for the particular choice £ — a + b, therefore c — and 2£ — 2a — b < k, i.e., b < k, 

Xpo,i3 = 0, \P\ < k. 

But that says that A ± n</;. 

In what follows, it is convenient to consider, for any sequence or indexed 'set' (vj : 
j e J) of vectors in some linear space Y over the scalar field F, the corresponding map 

with IFg denoting all scalar- valued functions on J with finite support. Note that the Vj 
enter the description of this linear map in exactly the manner in which the columns of an 
(m, n)-matrix A enter the description of the corresponding linear map IF"^ — > IF"^ : c i-^ Ac, 
hence it seems reasonable to call [vj : j e J] the column map with columns Vj. We can 
think of [vj : j £ J] as a row, much as we can think of a matrix as the row of its columns. 
The sequence [vj : j G J) is a basis for Y exactly when [vj : j G J] is invertible, in which 
case one might just as well refer to the latter as a basis for Y. Further, if (Ai : i G /) is 
an indexed 'set' in the dual, Y' , of Y, then it is convenient to consider the corresponding 
map 

[Xi-.iel]' :Y ^W' {XiV : i e I), 

calling it the row map with rows Aj (hence the use of the transpose sign) since the Aj 
enter the description of this linear map in exactly the same manner in which the rows of 
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the transpose (i.e., the columns of A) of an (m, n)-niatrix A enter the description of 
the corresponding Unear map IF"^ : c ^ A^c. We can think of [A^ : i G /]* as a 

column, much as we can think of a matrix as the column of its rows. With this, we are 
ready to think, as we may, of the composition of such a row map with such a column map 
as matrix, i.e., the Gram matrix or Gramian 

[Xi-.ie lY[vj :jeJ] = {XiVj -. i e I,j e J). 

This use of the superscript ^ also seems consistent with the notation x^y := x{j)y{j) 

for the scalar product of x,y e JR'^, given that it is standard to think of the elements of 
JR'^ as columns. For completeness, we note that the composition [vj : j E J][Xi : i & /]* 
makes sense only when / = J, in which case it is a linear map from Y to Y and the most 
general such in case Y is finite-dimensional. 

Schaback [S] considers, for the n-dimensional space M spanned by evaluation at the 
elements of the given n-set X in IR*^, a basis A = [Ai, . . . , A„] graded (by degree) (he 
calls any such a 'discrete moment basis') in the sense that the sequence (k^ : i = 1, . . . ,n), 
with 

Ki := max{/c : Aj _L n<fc}, all z, 

is nondecreasing, and, for each k, [Aj : > /c] is a basis for M n _Ln<fc. 

One readily obtains such a basis from any particular basis [ni, . . . , Hn] for M, by 
applying Gauss elimination with row interchanges to the Gram matrix 

(9) (//,()" : ^ = 1, . . . , n, a e Ki) = [/XI, . . . , V, 

with the columns of 

V-.^ar-.aeTZi] [()--j = l,2,...] 

so ordered that j ^ \aj\ is nondecreasing. 

Indeed, Gauss elimination applied to an onto Gram matrix such as (9) can be inter- 
preted as providing an invertible matrix L (the product of a permutation matrix with a 
lower triangular matrix) and thereby the basis 

[Ai,...,A„] := [ill, . . . , iin]{L~'^y 

for M, and a subsequence (/?i, . . . , /3„) of {aj : j = 1,2,...) so that, for each i, the first 
nonzero entry in row i of [Ai, . . . , A^j^V^, i.e., in XiV, is the Pith. Since, by assumption, 
the map j ^ \aj\ is nondecreasing, this implies that 

= max{fc : Aj _L n<A:} = i = l,...,n, 

and, in particular, i Ki is nondecreasing. Further, if ^jC(i)Ai _L H^k, then, since 
Ki = \(3i\, also c{i)Xi _L n<fc. But this implies that 

[cii):\pi\<k]B = 0, 
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with the matrix 

B:= :\(5l\(5^\<k) 

square upper triangular with nonzero diagonal entries, hence invertible, and therefore 
c(i) — for \l3i\ < k. This proves that, for each /c, [Aj : > k] is a basis for M n -Ln<fc. 
Schaback then considers the polynomials 

Wj : X i-^ Xj \\x — •p'*-' , j = 1, . . . ,n. 

Note that, by (8)Lemma, 

deg Wj — 2Kj — K,j — Kj. 
This implies that the Gram matrix 

A^W = {XiWj : i,j = l,...,n) 
is block upper triangular since 

(AW)(i, j) - X,w, - (A, ® A,)||a: - yf^^ = (-1)'^^ (A,, A,-)«. 
is zero as soon as Kj > Kj. Further, for each k and with 

Ik := {i: Ki = k}, 

the diagonal block 

A'Wih, h) = (-1)^(A„ Xj)k : ^,i G h) 

is invertible, by (7) Corollary and the fact that [Xi : i E Ik] is a basis for an algebraic 
complement of -Ln<fe in _Ln<fe. We will use later that this implies that is a graded 
basis for 

F := ran := { a{j)wj : a e IR""} 
j 

in the sense that j i— > degtu^ is nondecreasing and, for each k, [wj : degWj < A;] is a basis 
for F n n<fc. 

For the moment, we only use the conclusion (which Schaback draws for the case that 
M is spanned by point evaluations) that A^W is invertible, hence 

P = P^ .- W{A^W)-^A^ 

is the linear projector that associates with each p e 11 the unique element f E F that 
matches p at M in the sense that 

/// = HP, lie M. 

Schaback also observes that Ps is of minimal degree in the sense that F minimizes 

degG := max{deg(7 : g E G} 

among all polynomial subspaces G that are correct for M in the sense that, for every 
pen, they contain a unique match at M. 

However, much more is true. Pg is of minimal degree in the strong sense (of, e.g., 
[BR2]) that it is degree- reducing, meaning that 

deg Psp < degp, p eH. 

This condition is shown, in [BR2] , to be equivalent to the following property more explicitly 
associated with the words 'minimal degree': 
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Definition. The Gnite-rank linear projector P onH is of minimal degree := 

dim(Gnn<fc) < dim(ranPnn<fc), A; e IN, 

for all linear subspaces G of 11 that are correct for ± ker P. 
Proposition. Ps is of minimal degree. 

Proof: Let G be a linear subspace of 11 correct for M. Then, G is n-dimensional 
and, for any bases A of M and W of G, respectively, the Gramian A^W is invertible. 

Choose, in particular, A to be a graded basis for M and to be a graded basis for G. 
Then, for any k G IN, the first dim(Gnn<fc) = ^{j : degWj < k} columns of the Gramian 
A^W = (XiWj) have nonzero entries only in the first ^{i : Ki < k} = n — ^{i : Ki > k} = 
n — dim(M fl _Ln<A;) rows. The invertibility of the Gramian therefore implies that 

dim(G n n<fc) <n- dim(M n ±n<fe). 

On the other hand, there is equality here when G — ranPg since, as observed earlier, 
Schaback's Wj form a graded basis while Ki = degWi, all i, hence n — dim(M fl ±n</s) = 
i^{i:Ki<k} = i^{j:degWj <k}. □ 

The proof shows that a linear projector P on 11 is of minimal degree if and only if 
dim(ran P fl Il<:k) + dim(± ker P ±n<fc) = dim ran P, all k. 

The polynomial interpolant Psp to p at M is, in general (see below), not the least 
interpolant PbrP of [BR2] to p at M, hence we are free to give it a name, and Schaback 
interpolant seems entirely appropriate (hence the suffix 5"). 

We now compare Pg and Pbr in the specific context of [S], i.e., when M is spanned by 
evaluation on some n-set X in IR^. In addition to being of minimal degree, each interpolant 
is constant in any direction perpendicular to the affine hull 

bX 

of X (or, fiat spanned by X), i.e., both satisfy 

Pfix) = {Pf){Pxx), 

with Px the orthoprojector of IR*^ onto \?X. 

For the Schaback interpolant, this follows from the fact that, for any x, and any 

yebx, 

\\x - y\f = \\Pxx - y\f + \\x - Pxx\f, 

hence, since Xj ± ^<Kj, 

Wj{x) = Xj{\\Pxx - -f + \\x - Pxxf)"' = ^j\\Pxx - -f "^ 

by (8)Lemma. 
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This readily implies that both interpolants coincide in case X is contained in a 1- 
dimensional Bat. 

Also, both projectors commute with translation, i.e., 

Ppi- + y) = iPp){- + y), 

and interact with any unitary change of variables as follows: 

(10) Pp{A-) = {Pp){A'-) 

for all real unitary matrices A, as follows for Schaback's projector directly from the obser- 
vation that, for all such A and any x and y, 

\\x-Ayf=\\AiA'x-y)\\'=\\A'x-y\\\ 

However, while (10) holds for P = Pbr and arbitrary invertiblc A (see, e.g., [BRl]), this 
is not in general so for P — Ps, due to the fact that it is the 'kernel' {x,y) i— > exp{x^y) 
on which Pbr is based rather than Schaback's {x,y) ^— > — 2/|p^). We would therefore 
expect the two interpolants in general to differ when dim bX > 1 and ran P isn't just some 

n<fe. 

The simplest such example occurs when X = {xi,...,X4}isa 4-set in IR^ that spans 
H^. In that case, the range of each projector is of the form 

n<2 + wJR 

for some homogeneous quadratic polynomial w. For the least interpolant, [BRl] gives w 
as 

X 1-^ X{x^ ■ )^, 

with 

4 

such that 

(11) ^a(j)a;, = 0, 5;]a(j) = 0. 

j 

Since X spans IR^ this says that A _L n<2 and, since ^X = 4, A is, up to scalar multiples, 
the unique such element of M. That means that Schaback's w is the leading term of 

X I— > A||a; — -11^. 

Since A _L n<2, and ||a; — = (II^IP ~ 2x^Xj + H^^^jP)^, this leading term is 

a;^ A(2||a;|p|| • f + 4{x^-f), 
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and this is, offhand, not just a scalar multiple of the least's w. For, there is no reason to 
believe that A|| ■ f = 0. 

E.g., with the specific choice 

X = (0,ii,i2,2;) 
involving ii := (1, 0) and '12 := (0, 1), we have 

A:/^/(^)-^(l)/(ii)-^(2)/(i2), 

hence 

All • f = Wzf - z{l) - z{2) = z{l){z{l) - 1) + z{2){z{2) - 1), 

and this is zero only for special choices of z. To be sure, it is zero when 2; = ii + i2, i.e., 
in case of gridded data, giving us then bilinear interpolation. 

To be sure, graded bases of the space M, of linear functionals at which to match given 
values by polynomials, have been used before in multivariate polynomial interpolation. For 
example, the multivariate 'finite differences' introduced and used in [SX] are easily seen 
to form such a graded basis. 

In particular, the construction of the least interpolant makes use of a graded basis A for 
M (constructed by a more stable variant of Gauss elimination, namely Gauss elimination 
'by segments') but obtains the polynomial space G of interpolants as the span of the 'leasts' 
of the Aj. To recall (e.g., from [BR2]), any A e H' is uniquely representable, with respect 
to the bilinear form 

IRH®n^IR:(/,p)^ D'^{f){0)Dy{0)/a\, 
by the formal power series 

CXD 

A:=J;aW, 

fc=0 

with 

AW:= ^(A(n(r/a!, all /c. 

|a|=fe 

Then 

K := max{fc : A _L n<fc} = mm{k : X^''^ ^ 0} 

is known as the order of A, and At**' is, by definition, the least of A (with k taken to be 
— 1 when A = and, correspondingly, A^"^] — 0). It is easy to see that G = span{A''**' : 
i = 1, . . . ,n} depends only on M and not on the particular graded basis A for M used, 
and is spanned by homogeneous polynomials, hence, equivalently, is dilation-invariant. For 
the most striking properties of G (such as a finite list of constant coefficient differential 
operators whose joint kernel is G), see [BRl] or [BR2]. 

Acknowledgements Thanks are due to Tomas Sauer for a constructive reading of 
what I thought was the final draft. Further, on receiving a preprint of the present note, 
Robert Schaback informed me that, in the meantime, he, too, had extended (5)-(8) to 
arbitrary linear functionals, albeit with different proofs. 
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